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Abstract 

We consider weak measurements (WM) of a spin observable in quantum mechanics beyond the 
usual asymptotic regime. This is done by obtaining the exact time-dependent wave functions of the 
measuring apparatus for general non-ideal measurements. Ideal strong measurements and the usual 
WM regime are obtained as the two extreme limiting cases of our exact treatment. We show that 
Q_| in the intermediate regime non-ideal measurements lead to "semiweak" eccentric values that differ 

from the usual weak values. We further show that even in the WM regime the exact treatment 
leads to a meter behavior that can be markedly distinct from the one predicted by the usual WM 
formalism. We give several illustrations and discuss an application to the distinction of different 
realizations of the same density matrix. 

PACS numbers: 03.65.Ta 
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I. INTRODUCTION 



In recent times the issue of 'weak measurement '(WM) in quantum mechanics (QM) has 
gained a significant and wide interest in reahzing apparently counterintuitive quantum ef- 
fects. This path-breaking idea was originally proposed by Aharonov, Albert and Vaid- 
man(AAV) l| that introduced a measurement scenario in QM so that the empirically mea- 
sured value (coined as 'weak value') of an observable can be seemingly weird in that it yields 
results going beyond the eigenvalue spectrum of the measured observable. Since then, this 
idea has been enriched by a number of theoretical works ^-7\ and experimental works 8- 16| 
This novel technique of obtaining strange property of an observable has several irnplications, 
for example, to provide insights into conceptual quantum paradoxes |l8|^enti- 



fying a tiny spin Hall effect 



detecting very smal 



transverse beam deflections |l3| and 



19|. 



measuring average quantum trajectories for photons 

As it is well known a standard quantum mechanical measurement requires a strong cou- 
pling in order to produce an ideal one-to-one correspondence between the measured system 
and the device state (for this reason standard measurements are often coined as being strong 
or ideal). In contrast the WM scheme, as the name suggests, assumes a very weak coupling 
(so that the system state is kept grossly undisturbed) and interfering device states (so that 
the coherence in the device remains intact). To observe an effect of this weak perturbation in 
the device pointer's state, it is necessary to suitably select a particular subensemble which is 
technically termed as "post-selection". Given the initial state (known as "pre-selected" state 
in AAV's terminology) any 'strange' weak value of an observable can be observed if the post- 
selection is appropriately chosen. The weak values have another unusual property in that 
they can be complex. If the weak coupling approximation is made in position coordinate, 
such as, in a Stern-Gerlach (SG) setup, the real and imaginary parts of the complex weak 
value correspond to the momentum and position mean values of the post-selected pointer 
state respectively j^. In this case, if the weak value is real the effect can be observed in the 
momentum space distribution of the pointer state, i.e., the peak of the momentum distri- 
bution is shifted by an amount proportional to the weak value, possibly several times the 
corresponding eigenvalue of the observable in question. 

In the original AAV treatment, the weak measurement scenario is obtained as an approx- 
imate limiting case for weak couplings and nearly overlapping device states, whereas strong 
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projective measurements pertain to the opposite limiting case. However, between these two 
hmiting cases there is a whole continuous range of situations encompassing different regimes 



that have been scarcely investigated. Recently however, a few works |22N25| have studied the 
issue of weak measurement beyond the usual approximations made in the AAV formalism. 
The specificity of our work is that we analyze the general case of non-ideal measurements 
(followed by a post-selection) by starting from the exact solutions of the Schrodinger equa- 
tions coupling the system and the measuring device. The standard WM regime then appears 
for weak couplings and nearly overlapping meter states. But the other "non-ideal" measure- 
ment situations also display potentially interesting effects. We will show in particular that 
in these regimes not only can eccentric outcomes be obtained, but that the resulting "semi- 
weak" or "generalized weak"values can be tuned at will, allowing the device's pointer to be 
shifted more than the weak value. We will also see that it is possible to observe dichotomic 
outcomes of the pointer, as in the case of strong measurements, but with eccentric shifts, 
including in situations in which standard weak values cannot be defined. We will further 
show that even in the standard WM regime, some effects neglected in the AAV treatment 
can lead to exact weak values markedly different from the usual WM ones. All this features 
will be analyzed in a simple system that is analytically tractable - a spin- 1/2 particle passing 
through a series of Stern-Gerlach setups. 

The paper is organized as follows. In Sec. II we briefly recapitulate the essence of the 
standard weak measurement. In Sec. Ill we will introduce a particle with spin 1/2 passing 
through a series of Stern-Gerlach setups. The spatial wave packet will be considered as 
the probe that measures the particle's spin state. Solving rigorously the coupled solutions 
of the system allows to define naturally the non-ideal situation; the two limiting cases of 
strong and weak measurements will be explicitly given; the peculiar features visible in the 
intermediate regime - the semi-weak and the generalized weak values - will be described. 
Our findings are illustrated in Sec. IV, and an application to a quantum information task 
will also be considered. Our conclusions will be given in Sec. V. 

II. STANDARD WEAK MEASUREMENT SCENARIO 

The entire process of the AAV weak measurement procedure consists of three 

different steps: state preparation (usually termed as pre-selection), a strong projective mea- 
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surement for selecting a specific subensemble known as post-selection and in between the 
pre- and post-selection a weak interaction is introduced so that system state remains virtu- 
ally unaffected by this intermediate interaction (see Fig. 1). AAV demonstrated the entire 
process in terms of a series of three Stern-Gerlach (SG) setups. Let a beam of spin 1/2 
neutral particles, say neutrons, pass through the SG setups. The first SG is used to prepare 
the spin pre-selected state labeled by |xm)- The total initial wave function after the first 
SG setup is "^in = ipo{x)\xin) ■ The spatial part ipoi^) is taken to be a Gaussian wave packet 
peaked at the entry point (x = 0) of the second SG at t = 



(1) 



For simplicity the spatial part is written as being one dimensional - the spatial wavefunction 
along the y and z directions are trivial in that they are not affected by the Stern-Gerlach 
interaction given below. The initial momentum space wave function (poiPx) corresponding 
to Eq.p is 

The neutrons having the state = 'ipo{x)\Xin) then pass through the second SG setup 
that is used for measuring a spin observable, say, a^- The interaction Hamiltonian is given 
by if = f{t)fia.B where B = {bx,0,0) and fi is the magnetic moment of neutron. f{t) is 
a smooth function of t vanishing outside the interval < t < r and obeying f{t)dt = r, 
where r is the transit time during which the neutrons interact with the magnetic field. The 
total state after the interaction can then be written as 

*' = e —M^Mn). (3) 

The spatial part of the wavefunction can be considered in the context of the SG setup as 
describing the state of the measurement device: in a SG, the spin state is inferred from the 
wavepacket's deviation, that leads the particle to be found int the upper or lower planes 
when a strong measurement is performed. 

Here instead the magnetic field and transit time are taken to be very small, so that the 
system state remains essentially unaltered by the interaction. Then the exponential in Eq. 
([3]) can be expanded to first order. This is done in the original AAV treatment by performing, 
after the weak interaction, a strong projective measurement by using the third SG setup. 
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This "post-selects" the neutrons in a definite final spin state \xf)- The post selection allows 
to write the device state as 

'^fi.x) = (X/|e"^^^^'0o(a;)|Xm) = - + o(2) - ...|Xm)V^o(a:) (4) 



Neglecting the higher order terms leads to 

ijihrx {xMx\Xin) 



i^fi^) = iXflXin) 

which can be written as 



Mx) (5) 



where 

_ {XfWxlXin) 



^ iXflXin) 

i^fi^) = {Xf\x^n) e-*'^(^^)-^o(a^) (6) 

(7) 



(X/l Xin) 

is known as weak value of the observable ax- 

This weak value appears as a phase-shift in the pointer state's configuration space wave- 
function determines the position of the pointer in momentum space. Taking the Fourier 
transform of Eq. the pointer state in momentum space can be written as 

where p'^ = fibr. The final momentum distribution can then be written as 



(8) 



\MP^)\'=\{Xf\X^n)\'\MP^.-P>^U\' (9) 

where | (x/IXm) P is the probability of successful post-selection. Hence the final pointer 
position in momentum space is shifted by an amount {o'x)wP'x in contrast to the strong 
measurement case where the final pointer positions are shifted by ±p'^, corresponding to the 
eigenvalues ±1 of ax- 

From Eq.(I7]) it can be seen that {ax)w can be exceedingly large if pre-selected and post- 
selected state are nearly orthogonal. Rescaling the momentum to Px/p'x^ we see the pointer 
displays a broad distribution centered on {ax)w (or the real part of (cr^.)^ if the latter is 
complex), hence possibly well beyond the ranges of the eigenvalues. For instance for a 
specific pre-selected state, say, \xin) = \ te) = cos || tz) + sin || J,^) and post-selected state, 
say, \Xf) = |t)z) the weak value of ax is given by {ax)w = tan |. According to this formalism, 
{o'x)w can become arbitrarily large as ^ — tt, but simultaneously the probability of this post- 
selection, I {xflXin) P = litz I te)P = cos^ |, becomes vanishingly small (we will see below 



Weak 



Strong 



Figure 1: A series of Stern-Gerlach setup for implementing weak measurements of the spin operator 
Ox- The three SG setups account respectively for the state preparation, the weak measurement, 
and the post selection . 



from the exact treatment that both of these properties are erroneous, as the Umit ^ — )■ tt is 
well-defined). The final momentum distribution ([9]) after post-selection can be written as 



cos 



2\ 1/2 



vr 



r9 



exp 



^2 



(10) 



Note that in deriving the above relations we neglected the higher order terms in Eq. 
which is justified if for n > 2, 



(X/l(0-x)"|Xm) I << I (X/l 



{Xf\{^xT\Xin) I << {xMx\Xin) 



Moreover in order to resum Eq.([5]) into Eq.Q we need fibTx{ax)w/fi << 1; lastly noting 
that X is effectively governed by the spread 6 we may write the condition for the derivation 

as 

sf^iaX « 1- (11) 
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III. BEYOND THE STANDARD WEAK MEASUREMENT REGIME 



A. General remarks 

The usual weak measurement scheme that was presented in Sec. II involves rough approx- 
imations, some of which can usually be controlled (coupling value, width of device states), 
others, like the loss of unitarity due to the asymptotic expansion have consequences that are 
more delicate to tackle with in definite situations. Moreover there remains a whole continu- 
ous range of regimes for which the approximations made in the standard WM derivation do 
not hold, though the situation is still far from the ideal strong measurement limit. These 
intermediate non-ideal regimes that can be referred as 'semi- weak' situations are considered 
in the present Section. 

We first give the exact solution for a non-ideal Stern-Gerlach, that is solving the coupled 
Schrodinger equations for a wave packet of a spin 1/2 particle without making any specific 
assumptions for the coupling strength and the device states. We will then see how semiweak 
outcomes are obtained when a strong measurement is added after the neutron has emitted 
the non-ideal Stern-Gerlach. A crucial quantity that we introduce is the overlap / of the wave 
packets (recall the wave packets play here the role of the device states) that quantifies the 
strength of a given measurement situation. We will indeed see that / — yields the strong 
ideal measurement regime, while / — 1 corresponds to the usual AAV weak measurement 
scheme. For an arbitrary value of / the superposition of the device states results in a 
a distribution yielding the semiweak values. The different types of behavior of the pointer 
that can be obtained (one or several semiweak values, eccentric semiweak values or semiweak 
values falling between the eigenvalues, exact weak values that take into account details of 
the interaction Hamiltonian) will also be discussed. The illustrations are deferred to Sec. IV. 

B. Neutrons passing through a non-ideal SG setup: exact solutions 

The setup is the same as the one discussed in Sec. II (see Fig. 1). We omit the state 
preparation procedure and consider the initial spin state as given. Let a beam of neutrons, 
passing through the second SG setup be represented by the total wave function 




(12) 
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where = cos|| t^) + sin|| l^) is the initial state of the system (i,e., the spin). The 
spatial wave function ipo{'x) corresponds to a Gaussian wave packet which is initially peaked 
at X = at t = 0, given by 

where S is the initial width of the wave packet. The wave packet moves along the +y axis 
with the initial momentum Py (see Fig. 1). The inhomogeneous magnetic field^ B = {bx, 0, 0) 



27| is directed along the x-axis and confined between y = and y = d. The interaction 
Hamiltonian is Hi = fia.B where as above is the magnetic moment of the neutron. As 
the wave packet propagates through the SG magnet, in addition to the +y axis motion, 
the particles gain momentum along ± x-axis due to the interaction of their spins with the 
field. The time evolved total wave function at r (transit time of the peak of the wave packet 
within the SG magnetic field region) after the interaction of spins with the SG magnetic 
field is given by 

vl/(x,r)=exp f-'J^\^(^,t = 0) 



a 



^^+,(x, r) ® It), + /3^-x(x, r) ® |;)^ (14) 



where the device states ip^r^ (x, r) and ^/^.^ (x, r) are the two components of the spinor 

ip = ^ which satisfies the Pauli equation and ct = (cos | + sin |) and /3 = 

(cos I — sin I) . Note that Eq. (fT^ is an entangled state between position and spin degrees 
of freedom. The reduced density matrix of the system in the x-basis representation can be 
written as 

Ps=\ . (15) 



a(3r /3' 



where / is the overlap 



1= [ r+A^,r)ij.A^,T)d^x. (16) 

J V 

that quantifies the weakness of the measurement. The inner product / is in general complex 
but here in our case I is always real and positive. The value of the I can range from to 1 
depending upon the choices of the relevant parameters, such as, the degree of magnetic field 
(6), the width of the initial wave packet (5) and the transit time through the field region 
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within SG setup (r). We calculate the analytical expressions of i/jj^^ (x, r) and i/j^^ (x, r) by 
solving relevant Schroedinger equations. 

The two- component Pauli equation for and ip^x can then be written 



ih 
ih 





h^ 


dt 


2m 




h^ 


dt 


2m 



(17) 
(18) 



The solutions of the above two equations at t = r upon exiting the SG are as follows (for a 



detailed derivation, see Ref. 



20|) 



3" 6Xp 

(27r52)4 



X exp 



2m ■ 



4(52 



PyT\ 


Py 


p'x^ 


2m J 


J- 


' ft IJ 



(19) 



x:r 



^exp 

(27r(52)4 



X exp 



-A 



_ IV^\ Py_Px 
2m J h 



h 



(20) 



where A 



the evolution. 



p'^ = fibr, and the spreading of the wave packet is neglected throughout 



Here ip^^ (x, r) and tp^x (x, r) representing the spatial wave functions at r correspond to 
the spin states |'|")^ and \\)^ respectively, with the average momenta (p)^ and {p}i, where 
{p)u ~ i^P'x^Py^^)- Within the magnetic field the neutrons gain the same magnitude of 
momentum p'^ = fibr but the directions are such that the particles with eigenstates | 
and I 4)x get the drift along +x-axis and —a;- axis respectively, while the y-axis momenta 
remain unchanged. 

From these analytical expressions of ip+x (x, r) and ip^^ (x, r) given by Eqs. (fT9l) and (l20|l 
it is straightforward to compute the inner product / [Eq. f|T6|) ] given by 



exp 



(21) 



8m2(52 fi^ 

which explicitly depends upon the choices of the parameters b, 6 and r. 

Now, after emerging from this non-ideal SG magnet the neutrons represented by the 
entangled state given by Eq. (ll4l) enter to another SG setup where a strong measurement to 
be performed and the neutrons are post-selected in a specific spin state. 



C. Subsequent strong measurement: post-selection and final pointer state 



For this purpose, we consider immediately after the wavepacket exits the WM Stern- 
Gerlach a subsequent strong measurement of the spin observable a^. The knowledge of the 
exact solutions allow to treat the strong and weak measurements on the same footing. The 
approximate magnetic field in this case is B = (0, 0, b'z) and T is the time during which the 
interaction occurs. A strong measurement requires the magnetic field b' to be sufficiently 
strong so that the spatial wavepackets emerging from this SG setup are orthogonal, i.e., the 
relevant inner product, analog to Eq. fl21| vanishes. If the neutrons having the position-spin 
entangled state given by Eq. (ll4l) enters to final SG setup then the time evolved state that 
exits from the setup can be written as 



*(x,r + T)= ^V'+x+.i 



a 



x, r + T) It), + ^V'+x-H^- 



a 



X, r 



(22) 



+ ^V-x+.(x, T + T) It), + ^ + It), 



The inner products {ip^x+zl '^+x-z) and {ip-^^zl "^-x-z) at t = r + T vanish for strong 
measurements. The states 'ip±x±z{^, t+T) can be calculated by using ?/'+^.(x, r) and ?/^_x.(x, r) 
given by Eqs. (fT9|) and (pO|) respectively as the initial position wave functions. Note here 
that since the inhomogeneous magnetic field is only along the 2;-axis, within this SG setup, 
the neutrons move freely along the two other directions. Hence the analytical forms of the 
states ip±x±z{'^, T + T), are given by 



ip+x ±z (x; T + T) 



1 

^exp 

(27r(52)3 



i^T^i^Y + {y 



2m J 



X exp 



-A - A' ± 



h 



452 

Vy{.T + T)\ py 



2m 



, Px 

h h 



X 



2m 



(23) 



ijj-^ ±z (x; T + T) 



1 

3- exp 

(27r(52)3 



i^T^ifr+{l 



X exp 



n 



452 



y _ Pyi'r + T) \ Py L I P'xT 
2m J h h \ 2m 



(24) 
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where A' = and p'^ = fib'T 

Assume we post-select the neutrons having spin state | f^) then the post-selected position 
space wave function is written using Eq. (p2l) : given by 

<fpost (x, r + T) = [a^+.+z (x, r + T) + (x, r + T)] (25) 

The states -i/'+^j+zlx, T + r) and '?/;_2._|_^(x, T + r) are separable in x,y,z and can thus be 
written with obvious notation as 'ip±x+z{.^, T + r) = ■ip'^^{x, T + T)il'y{y, T + r)?/'+z(2:, T + r). 
Eq. (l25l) becomes 

^Po,i(x, r + T)~ ^pyiy, T + r)V^+,(z, T + r) T + r) + T + r)] (26) 

We see the y and 2; dependent parts of the total wave functions do not play any significant 
role in this context are thus integrated out. The motion along the x— axis within the final 
SG setup is free so we can neglect the T dependence of iIj'_^^{x,T + r). Thus the final 
post-selected pointer wave function, depending only on x can be written as 

%osti^, ^) = ^ H'+xi^^ ^) + r)] (27) 



where 



i^±Ax,T) = ^^exp 

(27r(52)3 



452 h 



(28) 



obtained from Eqs.([23D and (El]). 

The corresponding momentum space wave function is 



^postiPx, ^) = [«0+x(Px, r) + /30_^(p^, r)] (29) 



where (j)±x are obtained by taking the Fourier transform of Eqs. fl28l) . yielding, with given 
by Eq. © 

..(p.,.)^.o(p..P>.p(-Sl.^-^) (30) 

The momentum space solution ( 129|) is general in that no restriction on the coupling constant 
or on the width of the probe states have been introduced. We can now look at the device's 
momentum distribution in different regimes that will be characterized by the value of /. 
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D. Strong measurement limit: 1^0 

Let us now assume we tune the relevant parameters {b, S or r) so that / — )■ 0. In this 
case the reduced density matrix given by Eq. fllSp of the system becomes 



which is diagonal implying the strong measurement limit in that there is a one-to-one cor- 
respondence between the pointer state (resp. and the spin state (resp. 
Hence, in this case we expect two spatially separated peaks pointing at the eigenvalues of 
the spin observables ax- In Fig. (2a) we depicted the device momentum distribution given 
by Eq. (l29l) where we find the peaks at ±p^. 

E. Standard weak measurement limit: 1^1 

The WM limit is just the opposite limiting situation of strong measurement. From Eq. fl2ip 
we can see for example that for a fixed 6 we can choose the other parameters b and r such 
that / ^ 1 is obtained. The reduced density matrix of the system becomes 



i.e., the coherence in the system is considered to be mostly unaffected. There is no correspon- 
dence between the system and the device states. This leads to the usual weak measurement 
scheme proposed by AAV. Indeed let us start with the pointer position space wave function 
given by Eq. (l28l) since the coupling is between the spin and the position variable. If the 
parameters b and r are sufficiently small the value of I could be close to unity for a fixed 6. 
In this case we can neglect the higher order terms of br keeping only first order term. We 
can then write the states ip'_^^ (x; r) given by Eqs. (l28!) as follows 




(31) 




(32) 



In this limit Eq.( l27I) becomes 



V2 



[a (1 + ip'^x/h) + 13(1- ip'xX/h)] 



(33) 
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Putting the values of a, /3 and ipo{x), and simplifying we get 

^, , . cos I / .p'xtan|\ , 

« -p (-55 + '^x^ j (34) 

The Fourier transform of Eq. flM|) gives the pointer wave function in momentum space yield- 
ing the distribution 

l^/(Px)r - cos - ( — 1 exp 
which exactly matches the Eq. ([TO|) of the standard WM formulation. 

F. Nonideal situations: semiweak and exact weak values 

When the overlap / lies in the intermediate range < / < 1 an imperfect coherence 
remains in the system state. After post-selecting on a given final spin state, the pointer 
distribution is obtained from Eq. (127|) . therefore displaying an interference between the 
overlapping meter wavefunctions correlated with orthogonal spin states. 

The pointer behavior depends on the details of the interference: by monitoring /, eccentric 
pointer shifts can be obtained in this intermediate range. We coin these resulting shifts as 
'semi-weak' values. Depending on the resulting interference, that can be tuned by changing 
the values of the parameters (such as b, 6 or r) the pointer distribution can display a single 
maximum, akin to the usual weak measurement regime with the additional feature that this 
maximum does not necessarily appear at the AAV weak value given by Eq. ([7]), but can be 
chosen by fine-tuning the parameters. Other behaviors of the pointer can be obtained, in 
particular profiles with two maxima, similar to the two peaks that characterize the strong 
limit, but with the maxima shifted far from the eigenvalues. 

Even in the usual / ^ 1 weak measurement regime pointer shifts can be obtained in cases 
in which the usual WM formalism does not predict any, such as in situations in which the 
pre- and post-selected states are orthogonal. More importantly when the oscillating terms in 
Eqs. f l28l) play a role (that is when the wavelength of the oscillation becomes of the order of 
magnitude of the distribution width), the meter distribution becomes significantly different 
than what is predicted by the usual WM formalism. The features produced by these 'exact' 
weak values will be illustrated below. 



2'5^(Px-p;tan|) 



(35) 



13 



It is interesting to note that several recent works |22h25 
the usual WM formalism in the weak regime. Wu and Li 
to second order (one order beyond AAV's derivation) and showed this was sufficient to treat 



have also attempted to go beyond 
23| took the asymptotic expansion 



the orthogonal pre and post-selected states case. These results were extended in Ref. |24| 
where it was shown that the asymptotic expansion can be resummed when weak values of 
projectors are considered. A different but equivalent approach 22|, |25| follows the early 
method of Sudarshan and collaborators 2I consisting in formally expanding the preselected 
state in terms of the eigenstate of the weakly measured operator. This allows to recover the 
standard WM in the limit of weak couplings. The main difference of all these approaches 
with our present work lies not only in the method employed - we explicitly solve the full 
Schroedinger equation of the entire system (measured subsystem and pointer apparatus) -, 
but also in the fact that by doing so, we take into account the motion of the system and of the 
pointer, generated by the respective self-Hamiltonians. The strong, weak or non-ideal nature 
of the measurement is then seen to depend on characteristics of the full solution, not only 
on the coupling strength. For the weak measurement of spin in a Stern-Gerlach developed 
here, the different parameters that determine the non-ideal character of the measurement 
were encapsulated in the quantity /. 



IV. ILLUSTRATIONS AND APPLICATIONS 



A. Semiweak value distributions 



We give here an illustration of the properties of the meter state distributions in the non- 
ideal regime giving rise to the semiweak values examined in Sec. III. For definiteness the 
examples given will correspond to the explicit formulas given above, involving the initial 
wavefunction Eq.([T]), with the initial spin in state \t)g, an intermediate weak or semiweak 
measurement of the system (spin) observable cr^., and a final post-selection along \t)z- 



1. Strong and weak limits 

Consider ffist the strong and weak limits. A typical example is given in Fig. 2, where 
the device state distributions are plotted in configuration (x) space as well as in momentum 
(P) space (recall the only non-trivial axis is the one along x, the WM axis). For the strong 
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Figure 2: (Color online)The momentum distribution | ^pn^^ (p-r -t) p (Eg I29p is plotted when: (a) 
/ ~ (ideal projective measurement situation); and (b)/ ~ 1 (weak measurement situation). The 
values of relevant parameters b = lOOG/cm, r=1.4xl0~^s, and 5 is 1cm and l/50cm for (a) and 
(b) respectively. The weak value is (o"a:)^ = 16.2 for 0=173.5°. 

case the post-selection along is applied after a first strong measurement along that 
couples perfectly each wavepacket (meter state) with the the relevant spin state or H)^. 
In X space this yields a single broad structure at the eigenvalue of the az measurement, 
whereas in P space two peaks are visible(see, Fig. 2(a)), each corresponding to the direction 
of the momentum along the a;-axis: the strong measurement condition correlates perfectly 
this momentum with the spin state. 

In the weak limit, a single peak is still visible in x space, but in P-space the two inter- 
fering Gaussians overlap almost perfectly (/ ^ 1). The result is itself (approximately ^) a 
Gaussian, as predicted by the WM formalism, but with a maximum displaced at the weak 
value p',j. tan ^ (see Fig. 2(b)). The meter has therefore a broad distribution centered at a 
value several times higher than the spin eigenvalue with the single peak pointing here at 
{a^)u, = 16.2 for 6 = 173.5°.. 

2. Semi-weak values 

We first illustrate in Fig. 3 the generic situation as / is varied from the ideal strong 
situation (/ ^ 0) to the weak limit (/ ^ 1). We see that the meter distribution goes 
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Figure 3: (Color online)The momentum distribution |^p^cf(ff^. r)|^(Eq l29p is plotted for four dif- 
ferent values of the inner product(I) depending upon the suitable choices of the parameters. Here 
(a) and (d) represents the strong and weak measurements, and (b)-(c) correspond the semiweak 
measurement situations. The weak value is {(Tx)w = 13.3 for 9 = 171°. 

smoothly from the strong type of behavior to the one that characterizes the weak limit. 
Note however that the maximum does not appear anymore at the eigenvalues nor at the 
usual weak value Eq. (jT}) but at the semiweak values, that is the maximal values of the 
distribution obtained from Eq. (l29l) . In Fig. 3(c) for example the left peak is shifted to 
— 14.2 while the weak value is {ax)w = 13.3 (shown in Fig. 3(d)). This semi-weak value 
can be chosen at will by tuning the different parameters; approximate analytical expressions 
may be obtained by expanding \^post{PxiT)f in terms of the relevant small parameters (r, 
b or a combination thereof). 

3. Exact weak values 

A different feature of the exact formalism developed above concerns the existence of 
a meter distribution peaked at eccentric values in the usual weak regime (as far as the 
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weakness of the coupling and the meter states are concerned), but for which the usual AAV 
formalism does not yield any result. For example when the initial and post-selected states 
are orthogonal the usual weak value given by Eq. ([7]) is undefined. This undefiniteness stems 
from partially expanding the exponential after the post-selection [see Eq. (^]. An exact 
treatment does not have this problem. Fig. 4 shows the exact meter momentum distribution 
for orthogonal pre and post-selected states, for parameters near the 1^1 limit. In this 
situation, the post-selected meter state in momentum space is given by 

^postiPx, ^) = ^ [(p+x{Px, r) - (t)-x{Px, t)] (36) 

The analytical form of the momentum distribution near the / ^ 1 limit is readily obtained 
(up to a constant factor) from the Fourier transform of Eq. (l33|) as 

\^post{Px,r)f cx plexp (37) 

This distribution has two peaks at px = ±h/\/26 as illustrated in Fig. 4. The type of 
feature illustrated in Fig. 4 was first obtained theoretically in an approximate form (by 
neglecting the evolution of the pointer inside the Stern-Gerlach) by Duck et a/.Q, and more 



recently revisited in some of the works 



25[ going beyond the standard WM (eg in Ref. 
23| a formula for orthogonal pre- and post-selected states is obtained by consideriiig the 



interaction up to second order). Experimentally the feature was observed quite early 8|. In 



a more recent work|26|. the electric field analogue of Eq. (l37|) was derived in the framework 
of a phase-amplification technique using classical wave optics. Note that the exact weak 
value is perfectly finite whereas from the definition given by Eq.([7]) one might be tempted 
to conclude (incorrectly as this would violate the conditions given by Eq.(8) under which 
the AAV approximation holds) that the weak value can be arbitrarily large as the pre and 
post-selected states become orthogonal. 

Another instance in which the 'exact' weak value differs from the one defined in the 
asymptotic treatment appears when the oscillating terms of the exact solutions (Eq i30|) 
play a role. This happens when the wavelength of these terms become comparable to the 
width of the momentum distribution. A first example is shown in Fig. 5, when the pre- 
and post-selected states are identical. Then according to its definition ([7]) the weak value 
should simply be given by the average of the weakly measured observable in the initial 
state. It can be seen in Fig. 5 that this rule is spoiled by the oscillating terms, that create a 
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Figure 4: (Color online)The momentum distribution (p^. r)p(Eq l36|) is plotted when / w 1 but the 

pre- and post-selected states are orthogonal - a situation in which the standard weak value is undefined. 




Figure 5: (Color online)The momentum distribution \^postiPx,T)\'^ [Eq. ([29| ] is plotted when the pre- 
and post-selected states are identical. The values of relevant parameters b = O.OOlG/cm, t=1.4x10^^s, and 
5 is 5 X 10^^ cm . 

multiple peak structure with eccentric maxima. By tuning the parameters entering the exact 
wavefunction, the wavelength of the oscillating terms can be reduced, yielding a momentum 
distribution with multiple peaks, as shown in Fig. 6. The origin of the multiple peaks is 
the interference between the two meter states emerging from the second SG setup along two 
opposite directions carrying the opposite phases. These phases (in momentum space) arise 
because of the separation of the two wavepackets in configuration space due to the spin- 
inhomogeneous magnetic field interaction occurring along opposite directions. That shows 
the necessity of the solution of the Schroedinger equation using the full Hamiltonian instead 
of taking only the interaction Hamiltonian or keeping only a few terms of the asymptotic 



expansion of the interaction coupling. Is 



the main interest in several works 



ote here that the average momentum that has been 



22h25| is not necessarily relevant when there are multiple 



peaks in the momentum distribution. 



18 




Pv 

(pj 

Figure 6: (Color online)The momentum distribution \^post{Px,T)\^ [Eq. ([29| ] is plotted when the pre- 
and post-selected states are identical. The values of relevant parameters b = O.OOlG/cm. t=1.4x10^^s, and 
S is 10^'^ cm . 

B. Application: Distinguishing between "identical" density matrices 

As an illustration of the practical usefulness of the exact scheme for weak values reported 
in this work, we give here an application to a quantum information task, namely the possibil- 
ity of distinguishing between "identical" density matrices. Consider the following situation: 
Alice prepares neutral spin 1/2 particles in some state, either or pz and sends them to 
Bob, whose goal is to guess the state. According to elementary quantum mechanics the spin 
density matrices 

P.^^\W{U + l\U){U (38) 

and 

p. = ^it.)(t.i + ^ii.)a.i (39) 

are identical and thus undistinguishable. We must therefore give an additional condi- 
tion: we assume Alice sends successive spins of alternate signs, ie Alice sends either 
^ = {\tx) , \ix) , It.;) , \ix) •••} or C = {It.) , It.) , It.) , \lz) •••}, each of the sets ^ and C giving 
rise to a specific realization of px or pz respectively. Bob must guess as fast as possible, that 
is by processing the lowest number of particles, whether Alice is sending ^ or (. 

With strong measurements, the best Bob can do is to measure either ax or and examine 
whether two successive measurements have identical signs. Suppose for instance Bob chooses 
to measure a^; the pointer displays two sharp peaks, one positive and one negative (see Fig. 
[7] (a)). Then if ^ is sent the sign of successive outcomes will alternate, whereas if ( is sent 
each outcome is equiprobable. So Bob's strategy will be to observe whether the consecutive 
measurements have alternating outcomes, in which case he will conclude it is likely Alice is 
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Figure 7: (Color online)Meter distributions (in arbitrary units) for the task described in Sec. IV. 2 
for strong [(a)], standard weak [(b)] and exact weak [(c)-(d)] measurements, with the magnetic field 
strength fixed in all cases at b = 0.02 G/cm and r = 0.07 s. In (a) ax is measured strongly {6=1 
cm) and the two peaks appear at the corresponding eigenvalues, irrespective of whether or is 
being sent, (b) shows the meter distribution corresponding to ^ (dashed, online blue curve) and 
to ^ (solid, online green curve) for a standard WM of dx (with 6 = 10~^ cm) and a post-selection 
along 55°. (c) shows the same situation displayed in (b) but with the exact WM solutions, (d) 
shows the configuration space meter distributions corresponding to the exact WM momentum space 
distributions shown in (c) with the same post-selection angle. When Alice launches ^ the left (online 
green) and right (online red) peaks are obtained for \lx) and |t^.) respectively, whereas for l^^) 
and 1^2) can both hit the left and right peaks. 

sending Indeed denoting by k the number of particles that have already been seen to be 
displaying an alternate series, the probability of continuing with this series for the [k + l)th 
particle if ( has been sent becomes 2~'' and decrease rapidly with k. 

With standard weak measurements the problem for Bob is that he ignores the preselected 
state: he must guess whether Alice is sending ^ or ^ by obtaining meter distributions that 
are very broad in momentum space, irrespective of the weak measurement and post-selection 
Bob chooses. Typically all the meter distributions have almost identical profiles, but different 
heights, which is how ^ and ( can be distinguished (see Fig. [71(b)). So despite the fact that 
the meter distributions will indeed be different, in practice Bob will need a great number of 
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particles in order to discriminate ^ from 

However by following the exact treatment in the non-ideal case, Bob can set (by changing 
the SG magnetic field strength and passage time) the interference between the meter states 
in momentum space so that the detection of ^ and ( result in totally different probability 
distributions. Moreover by weakening the coupling constant, Bob can still keep markedly 
different profiles in momentum space (see Fig. [7](c)) while obtaining non-overlapping profiles 
for the meter in configuration space. The advantage is obvious: one has exactly the same 
information as the one obtained with strong measurements (but in configuration space, see 
Fig. [7](d)), and in addition as more particles are detected the detection curve in momentum 
space unambiguously reveals whether ^ or is being sent. It is crucial to note that the two 
distinctly positioned peaks visible in configuration space are a feature of the exact solutions 
- in the standard WM formalism the configuration space wavefunctions are identical up 



to some global factor. In practice 



2l|, each odd numbered result is registered separately 



from the even numbered events, and this is done for each of the two post-selected states. 
The procedure can be stopped when enough events are registered so that the curves can be 
discriminated by employing standard statistical tests minimizing the distance between the 



expected curve and the obtained outcome 



2l|. 



V. SUMMARY AND CONCLUSIONS 



In this work we have derived the exact behavior of a measurement apparatus for arbi- 
trary non-ideal measurements. This was done in the analytically tractable case of a spin 
measurement in Stern-Gerlach setups, for which the exact time-dependent wavefunctions 
can be obtained. This has allowed us to investigate the validity of the usual AAV weak 
measurement formalism, both in the intermediate regime, far from the usual WM limit, 
but also in the WM limit when the specific features of the system-meter interaction are 
explicitly accounted for. In doing so we have been led to introduce 'semiweak' and 'exact 
weak' values that describe the meter distribution and the associated eccentric values in these 
circumstances (including, but not limited to, the case of orthogonal pre- and post-selected 
states, for which the usual weak- value is undefined). 

We have also seen through our exact approach that including the full Hamiltonian (rather 
than restricting the treatment to the interaction Hamiltonian) may have crucial consequences 
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on the behavior of the eccentric values. In this sense, an exact approach brings additional 



features than those that can be retrieved by other recent works |22|-|25| dealing with ex- 
tending the original weak measurement approach: it is valid irrespective of the interaction 
strength and it yields fine details due to the specificities of the interaction. On the other 
hand, the practical applicability of an exact approach hinges on the ability to solve (ana- 
lytically or numerically) the full Schrodinger equation, a task that can become difficult for 
an arbitrary system. In all cases, going beyond the usual WM formalism as well as the rel- 
evance of performing exact calculations for a given experimental situation should be helpful 
in analyzing experimental results obtained in nonideal measurements. 
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20|. However on average 



We need at least another component to make it divergence free 
the effect of these additional field components can be neglected under proper circumstances, 
resulting in this effective field usually found in textbooks and also employed in all the previous 
works dealing with weak measurements of spin. 
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